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Abstract 

The symmetry conditions for the devefopment of a macroscopic electrical polarization as a sec- 
ondary order parameter to a magnetic ordering transition, and the constraints on the direction 
of the polarization vector, are determined by a non-conventional application of the theory of ir- 
reducible co-representations. In our approach, which is suitable for both magnetic and structural 
modulations, anti-unitary operators are employed to describe symmetry operations that exchange 
the propagation vector k with -k, rather than operations combined with time-reversal as in clas- 
sical corep analysis. Unlike the conventional irreducible representations, co-representations can 
capture the full symmetry properties of the system even if the propagation vector is in the interior 
of the Brillouin zone. It is shown that ferroelectricity can develop even for a completely collinear 
structure, and that helical and cycloidal magnetic structures are not always polar. In some cases, 
symmetry allows the development of polarization parallel to the magnetic propagation vector. Our 
analysis also highlights the unique importance of magnetic commensurability, enabling one to de- 
rive the different symmetry properties of equivalent commensurate and incommensurate phases 
even for a completely generic propagation vector. 



PACS numbers: 75.25.+Z, 77.80.-e, 61.50.Ah, 73.22.Gk 



I. INTRODUCTION 



There has been a recent surge of interest for improper ferroelectric transition-metal com- 
pounds where the onset of electrical polarization is inducved by a transition to a complex 
magnetic state-i. Much of the discussion in the literature has focused on establishing the 
microscopic mechanism that couples the magnetic moments with lattice displacements: sym- 
metric superexchange (sometime referred to as "super exchange striction") and antisymmet- 
ric exchange (the so-called Dzyaloshinskii-Moriya contribution) have both been discussed, 
often in the context of the same materials^i^, and shown in many cases to describe the details 
of the coupled magneto-electric transitions, including the change in sign or direction of the 
polarization in an applied magnetic field. The study of the crystal and magnetic symme- 
try of these systems represents an integral part of this work. Establishing the symmetry 
constraints upon the polarization vector for a given magnetic structure, regardless of the 
microscopic mechanism of magneto-elastic coupling, is extremely valuable. For example, if it 
is established that for a given magnetic structure the polarization vector P lies by symmetry 
along a particular crystallographic direction, the experimental determination of the direction 
of P can be used to corroborate or falsify the magnetic structure model. Conversely, the 
exact direction of P must be predictable from the microscopic coupling mechanism when 
symmetry allows P to lie in a plane or in a general direction. More generally, and in analogy 
with the well-know case of magnetostriction, if the broken symmetry allows the development 
of a macroscopic polar vector, one would always expect to observe ferroelectricity, provided 
that the measurement has sufficiently high sensitivity. The "new" class of multiferroic ma- 
terials we are discussing is in fact characterized by very small values of |P|. However, the 
magnitude of P must be predictable from microscopic models based on the specific observed 
spin arrangement, out of the many that are usually consistent with each symmetry class. 

A number of techniques exist to predict possible symmetries derived from a given parent 
structure through a given order parameter - a subject that was thoroughly developed in the 
1980's^. When the propagation vector of the structural or magnetic modulation lies in a high- 
symmetry point of the Brillouin zone, one can use standard irreducible representations^i^i^ 
(irreps) to generate the so-called image of the high-symmetry group - the finite set of 
matrices representing the group elements in the order parameter space for a given irrep. By 
analyzing the image one can determine the form of the Landau free energy as a function of 
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both primary and secondary order parameters, enumerate the possible invariance groups of 
the low-symmetry phase and establish compliance with the so-called Landau criterion for 
continuous phase transitions^. When dealing with magnetic transitions, some of the early 
work employed irreducible co-representations {coreps) instead of irreps. In this approach, 
first discussed by Wigner^, the time reversal operator in the black-and-white symmetry 
groups is antiunitary, as applied to the Schrodinger wavefunction describing the magnetic 
ground state. However, it was later noted (see for example^^) that the use of antiunitary 
operators is not actually necessary when dealing with classical spin systems, for which time 
reversal is represented by a simple sign change, and the use of coreps has been all but 
abandoned as a consequence. 

When the propagation vector of a magnetic or modulated crystal structure is not a high- 
symmetry point of the Brillouin zone, the situation is considerably more complex. First 
of all, the "mathematical" or "physical" representations are in this case always obtained 
by combining at least two little-group- zrrej>s with opposite propagation vectors (-|-k and 
-k). Secondly, the image of this representation can have an infinite number of elements if 
the propagation vector is incommensurate with the crystal lattice. Thirdly and perhaps 
more importantly, the symmetry properties of the modulated structure upon application 
of operators that exchange -|-k with -k, such as the inversion, do not emerge clearly, since 
these operators do not have an image matrix in either of the irreps. For these reasons, image 
analysis is not ordinarily applied to these problems. Instead, one usually resort to construct- 
ing Landau free energies containing mixed terms in the -|-k and -k order parameters, and 
examining the symmetry of the solutions. Harris and collaborators have extensively em- 
ployed this approach to study in some detail the specific case of incommensurate magnetic 
multiferroicsiiii^ii^i^i^i^. In particular, they have specifically pointed out the importance of 
inversion symmetry, the globally invariant forms of the Landau free energy and the situations 
where this symmetry is spontaneously broken, giving rise to a net polarization. 

In this paper, we follow a different approach to the problem of determining the point- 
group symmetry of a modulated magnetic structure and of defining the symmetry conditions 
imposed on the development of a macroscopic polar vector upon magnetic ordering with 
an arbitrary propagation vector k. In our analysis, we employ the mathematical tool of 
irreducible co-representations, but its significance is radically different from that of the 
"standard" corep analysis of magnetic structures. This difference is apparent when one 
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considers the role of antiunitary operators, that are conventionally employed to represent 
"black" elements of the groups inverting the direction of time. Here, we use antiunitary 
elements to describe operators that exchange +k with -k, as explained below. This approach 
enables one to employ a single corep and propagation vector instead of two irreps. Time 
reversal is not an essential ingredient of this method, which can be equally well employed 
to study incommensurate structural modulations (for magnetic structures, we treat time 
reversal classically, es explained below). Crucially, all the symmetry properties emerge 
naturally from this analysis, since the "little group" of the propagation vector is extended 
to include the inversion and all other operators exchanging +k with -k. This approach 
is not particularly new - it is implicit in the treatment of coreps described in the classic 
book by Kovalev- and its significance has been recently re-emphasized by Schweizer— . 
However, we do not believe that these techniques have been hitherto employed to determine 
crystallographic point groups and macroscopic observables, as we do herein. With respect to 
the approach followed by Harris and coworkers, the main difference is that we perform our 
analysis directly on the images, as in the case of high-symmetry k vectors. Knowledge of the 
Landau free energy form is therefore not required, making this method easier to implement in 
an automated and tabulated form, as appropriate for non-specialists. We nonetheless stress 
that the Landau analysis provides much more information than the point group of the low- 
symmetry phase, and is therefore the tool of choice for specialist theoreticians. Furthermore, 
our method is completely general - inversion symmetry is treated on an equal footing with 
other operators exchanging -|-k with -k. Our only restriction is that the magnetic structure 
be described by a single k, although the extension of our analysis to multi-k structures is 
quite straightforward. 

We show that ferroelectricity can develop even when the magnetic structure is described 
by a single order parameter, and that P || k is allowed by symmetry in some cases. Fur- 
thermore, our analysis evidences the crucial difference between incommensurate and com- 
mensurate magnetic structures even for propagation vectors inside the Brillouin zone, and, 
for the latter, shows that the global phase has an influence on symmetry. This is partic- 
ularly counterintuitive, since the global phase affects neither the magnetic energy nor the 
intensity of the magnetic Bragg peaks in (unpolarized) neutron diffraction. In fact, the ob- 
servation of a non-zero electrical polarization can be used to discriminate between otherwise 
indistinguishable magnetic structures. 
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The paper is organized as follows: in Section [TTl we describe the use of corep to determine 
the point-group symmetry of a magnetically modulated system (the extension to lattice 
modulations is straightforward, and will be described elsewhere). In Section HIT] we describe 
in detail the application of this method to a number of topical multiferroic systems with 
commensurate or incommensurate magnetic structures. Section IIVI contains summary and 
discussion of the results. For completeness, in Appendix |A] we provide a brief overview of 
the theory of irreducible co-representations. 

II. THEORY 

As for all macroscopic observables that are even by time reversal, the existence of an 
electrical polarization and the restrictions on its direction are defined by the structural 
point group S of the magnetically ordered structure m. If we write a generic element of the 
paramagnetic space group G in the form g = {r|w} (Seitz notation), where r is a proper or 
improper rotation belonging to the paramagnetic point group and w is a translation, then 
r G S" if and only if there is an element g E G for which gm = ±m. If we only consider the 
representative elements of the group, go = {r|v}, where v is a non-Bravais translation and 
rotations appear only once in the representative set, then 

r G 5 < > {r|v}m = ±{E\t}m (1) 

where t is a Bravais translation and E is the identity. In other words, the corresponding 
representative element must be equivalent to ± a lattice translation. Here, we deal with 
the time reversal symmetry by considering both positive and negative eigenvalues, rather 
than combining time reversal and complex conjugation as it is often done^ii^. It can be 
shown that the two approaches are completely equivalents^. It is important at this point 
to recognize that the magnetic structure m is real, i.e., it is not a generic element of the 
linear space V defined over the complex field as a collection of axial vectors associated with 
atomic positions in the crystal. If m is described by a single propagation vector k, we can 
always write m = e^*'' *'?/' + e~*'^'*'?/'*, where e^^^''^ip is a generic element of the complex- valued 
subspace of V associated with the "arm" k. In the most general case, e^^'^'^ip and e"**'*?/'* 
transform with distinct (albeit complex-conjugate) representations of G, and in fact may 
not even belong to the same star. Conventional analysis using the irreps of the little group 
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Gk (i.e., of the group of operators leaving the propagation vector invariant) deals with the 
two Fourier components separately, combining them later in a single " physically irreducible" 
representation of higher dimension. As we shall see, this method is unable to capture the 
full symmetry properties of m. In our specific case, this means for example that m can 
be centrosymmetric even when the inversion operator / does not belong to the magnetic 
little co-group, as it is always the case if k is a non-Lifshits vector (2k ^ L* , L* being 
the reciprocal lattice). On the other hand, working with m (i.e., the two representations 
simultaneously) is extremely inconvenient, since, unlike e~^^^''^ip, m is not an eigenvector of 
the pure translations. With this in mind, it is useful to reformulate Eq. [1] as a condition on 
g+ik t^ rather than on m: 



r e S < > {r|v}e+ V = ± < 



g-ik-tog+ik-t^ 

or (2) 

g+ik-tog-ik-t^* 

Eq. [2] can be further simplified by introducing the operator of complex conjugation K: 



r G 5 < > < 



or (3) 

In other words, e~^^^'^il) must be an eigenvector of the operator {r|v} or of the operator 
ii'{r|v} (or of both), with eigenvalues corresponding to ± Bravais translations. It is worth 
pointing out already at this stage the crucial difference between incommensurate and com- 
mensurate propagation vectors. In the former case, any eigenvalue will do, since a suitable 
translation can always be found to equate any phase in the exponential. On the contrary, 
only a finite number of phases are available in the commensurate case. Based on Eq. |3l 
it becomes natural to consider the mapping not of the group G, but of the direct product 
group {E, K} C?) G. Here, the subtlety is that the images of elements of the form Kg must 
be antilinear and antiunitary operatora^S. it is noteworthy that the usefulness of this ap- 
proach is by no means limited to magnetic structures, and is equally applicable to structural 
modulations, provided that we consider only the " -|-" sign in Eq. [1] - [31 Homomorphisms 
of {E,K} ® G are known as co-representations (coreps) of G, and their theory has been 
extensively developed^^^^ (see Appendix |X] for a summary of this theory). In essence, the 
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corep analysis consists of 3 steps: 

1. Determination of the the subset of {E, K}®G that leaves k invariant, the equiva- 
lent of the "little group" G^. contains all the operators of the conventional "little 
group" , plus operators of the form Kg, where K is the complex conjugation and g E G 
exchanges -|-k with -k. Consequently KI G if J G G. 

2. Determination of the coreps and their image matrices. The complete analysis has been 
done by Kovale\*^ for all space groups and k vectors, and all is required is to refer to 
the tabulated values therein. 

3. Determination of the characteristic (basis) vectors for each corep. This can be done 
directly, by applying a projection method similar to the standard irreps, or, perhaps 
more easily, by symmetrizing the irrep basis vectors, as explained inil. One must 
keep in mind that, unlike the case of irrep basis vectors, corep vectors cannot be 
multiplied by an arbitrary complex constant, because of the antiunitary character of 
the associated operators (see below). 

Once this analysis is done, the symmetry condition in Eq. [3] can be thoroughly explored by 
determining the spectra of the unitary and antiunitary operators (images) associated with 
the various coreps. Crucially, operators such as the inversion / that exchange -|-k with -k 
(not included in the conventional irrep analysis) will now be represented by their antiunitary 
counterparts (e.g., KI), which do possess an image. Spectra and eigenvectors for the uni- 
tary operators are found in the usual way by diagonalizing the corresponding matrices. The 
method to determine the "characteristic vectors" of an antiunitary operator A is described 
by Wigner— . In particular, it is shown how to construct a full set of orthonormal vectors 
Vi . . .Vn that are invariant to both A and the unitary operator A"^, by linear combinations 
of the the eigenvectors of with eigenvalue = +1. Linear combinations of the Vj^s with 
real coefficients are also invariant by A. Multiplication of Vj (or of a real-coefficient linear 
combination thereof) by a phase factor e*"^ results in an eigenvector of A with eigenvalue 
A = e^^*'^. Linear combinations of eigenvectors with complex coefficients are generally not 
eigenvectors. A full spectral analysis of each operator is often not necessary, particularly 
when the aim is to establish the symmetry of an experimentally determined magnetic struc- 
ture (see examples below). 
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III. EXAMPLES 



In this section, we will analyze the symmetry properties of some magnetic improper 
ferroelectrics from the recent literature, using the co-representation approach we described in 
the previous section. In each case, we will determine the matrix representatives (images) for 
the relevant propagation vectors/co- representations and the associated basis vectors for the 
magnetic sites. We will also determine the point-group structural symmetry for particular 
ordering patterns. 



A. Multiferroic behavior in i^E'MnOs 

The space group is Pnma (no. 62 in the International Tables^; we adopt the standard 
setting, as opposed to the Phnm setting used in some papers), and the propagation vector 
is (/i,0,0), with |U incommensurate or commensurate but generally in the interior of the 
Brillouin zone. This propagation vector is labeled as k-j in Kovalevli. We will employ 
the standard International Tables setting rather than the "Old Kovalev" setting (both are 
reported in Ref.— ). The small irreps for this space group and propagation vectors are all 
one-dimensional, and their matrices (complex numbers in this case) are reported in Tab. [H 

TABLE I: Small irreps (A) and coreps (D) of space group Pnma for propagation vector = 
(/u, 0, 0). The symmetry operators are in the Kovalev notation and correspond to the International 
Tables symbols hi = 10,0,0; /i2 = 2(i, 0, 0) x, i, i; hs = 2(0,^,0) 0,y,0; h^ = 2(0,0, i) 1,0, z; 
/i25 = 1 0, 0, 0; /i26 = n{0, i, i) i, y, z; /127 = mx,l, z; /i28 = ax,y,l; e = e+*''^ 

hi /l2 /l27 h28 
Kh25 Kh26 Kh3 Khi 



Ai/Di 1 e 1 e 

A2/D2 1 e -1 -e 

A3/A3 1 -e 1 -e 

A4/D4 1 -e -1 e 



All the coreps correspond to " Case a" described in Appendix |X] ; in other words, each 
irrep generates a single corep. In addition, the coreps can be set in diagonal form, as 
explained in Appendix |X1 The matrices of the antiunitary operators are equal to those of 
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corresponding unitary operators, as shown in Tab. HI The symmetry properties of each 
corep or combination thereof is now clear by inspection of Tab. [H while remembering that 
the antiunitary operators complex-conjugate all mode coefficients. In particular: 

1. Linear combination of corep modes with purely real or purely imaginary coefficients 
are always centric or anti-centric and cannot support ferroelectricity. 

2. Multiplication of a centric or anti-centric mode by a phase factor e^^ is always equiv- 
alent to a translation for an incommensurate propagation vector but not necessarily 
so in the commensurate case. For incommensurate propagation vectors, this analysis 
confirms that ferroelectricity cannot arise from a single magnetic order parameter—. 
As we shall see in the remainder, however, this is only true in general for 1-dimensional 
co-representations. 

3. Linear combination of two corep modes with arbitrary complex coefficients in general 
violates all the antiunitary operators and those unitary operators with different ma- 
trices for the two coreps. In general this will lead to the polarization vector being 
allowed in a plane containing the propagation vector. 

4. Cycloidal structures are the most important case, because they correspond to the 
magnetic structures proposed in the literature for the ferroelectric phases. When two 
components are summed in quadrature, they do not always violate all the antiunitary 
operators. All linear combinations of this kind, with the associated structural point 
groups and allowed directions of the electrical polarization, are listed in Tab. HTl 

TABLE II: Structural point groups for cycloidal structures of general formula aD^ + ibDp [a and 
h are real coefficients). The allowed direction of P is indicated in parenthesis. A dot (.) means 
that the point group is non-polar, and no ferroelectric polarization can develop. 

Di D2 D3 D4 

iDi 2mm{x) 222(.) mm2{z) m2m{y) 

iD2 222(.) 2mm{x) m2m{y) mm2{z) 

iD^ mm2{z) m2m{y) 2mm{x) 222(.) 

iD/i m2m{y) mm2{z) 222(.) 2mm{x) 
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From Tab. [TTl we can see that the magnetic structure proposed by Kenzelmann^i, 
corresponding to the admixture D2 ^iD^, only allows a polarization in the y direction 
{z direction in the Pbnm setting proposed in Ref.Ji), as observed experimentally. 
The proposed magnetic structure is therefore consistent with the electrical properties, 
but no specific magneto-electric mechanism can be inferred from the observation. It is 
noteworthy that some combination of irreps induce non-centric, non-polar point groups 
(222 in this case). The absence of a center of symmetry should not therefore lead to 
the conclusion that ferroelectricity is allowed in some direction, nor the observation of 
a cycloidal structure lead to the conclusion that ferroelectricity is allowed. 

5. From this analysis, it is apparent that commensurate structures will in general have 
lower symmetry with respect to corresponding incommensurate ones - a well known 
general result^. Here, the obvious reason is that phase factors are not necessarily 
equivalent to translations. However, the symmetry may be higher for particular choices 
of the overall phase factor. An interesting example, which includes the magnetic 
structure proposed by Aliouane et al, described from the admixture: 

^ = e'"(aDi + ^6D3) (4) 

where a and b are real coefficients. This structure is always invariant by application 
of the mirror plane _L b (/i27)- Application of the two antiunitary operators i^/i4 and 
Kh26 yields (Tab. [I]): 

Kh26tlJ = Khi%l) = ee-'^'iaDi + ibD^) = e-'^^^+^f"^^ (5) 

Therefore the corresponding rotations belong to the structural point group only if the 
phase factor corresponds to ± a lattice translation, i.e. , if 

iv = -/Uvrn -|- 27rm (6) 

where n and m are arbitrary integers. If fi is incommensurate, Eq. [H] can always be 
satisfied to an arbitrary approximation. In this case, the structural point group is 
S = mm2, and the polarization is along the z axis. If /i is commensurate, only a 
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restricted number of phases are available, and Eq. [6] may be far from being satisfied. 
In this case, 5* = .m., and the polarization is in the x — z plane. The constant- 
moment "Aliouane" structure corresponds to this case with the particular choices 
fx = J, a = {2 + \/2)^, 6 = (2 — \/2)^ and cu = — |, for which Eq. [H] is satisfied. 
Consequently, S = mm2 and the polarization must be directed along the z axis. 
This is also the direction of the polarization found experimentally for the high-field 
phase of TbMnO^^. The microscopic model proposed by Aliouane, based on exchange 
striction, does in fact produce polarization, as required by symmetry. 

In i^E'MnOa, the Mn atoms are on centers of symmetry, and the application of the 
antiunitary operators does not generate more sites than those generated by the little group 
G^. Consequently, each instance of an irrep basis vector generates a single instance of 
the associated corep, spanning exactly the same subspace. Tab. IIIII lists the corep modes 
obtained by symmetrizing the conventional irrep modes, following the procedure described in 
Eq. IA3I It can be easily verified that the two modes generated by [A3] are linearly dependent 
via a single real coefficient. The magnetic structure proposed by Kenzelmann^^ corresponds 
to mx{D^) + imy{D2), which, as already remarked, only allows the polarization to be along 
the y-axis. A cycloidal structure of the same type but with spins in the x — z plane would be 
described as m^lD^) + imz{Di), which, according to Tab. [Til yields a polarization along the 
2;-direction (a;-direction in Phnm). This is consistent with the Ginzburg-Landau analysis 
performed by Mostovoy^^ for the specific case of cycloidal structures. It is important to 
remark that the direction of P as established by symmetry does not depend specifically 
on the direction of the magnetic moments or the type of magnetic structure (cycloidal, 
helical, etc.) For example, a helical structure of the type rUziD^) + imy{D2) and even some 
complex collinear structures (e.g., my{D-^) + imy{D2)) have exactly the same symmetry as 
the Kenzelmann^ structure. Naturally, magnetic measurements and neutron diffraction can 
all be used to distinguish between these possibilities, and to guide the analysis towards a 
microscopic model. 

The case of the RE sites is more interesting, because the atoms do not sit on a center of 
symmetry, and they are therefore split into orbits by the little group G^. The application 
of the antiunitary operators mixes the two orbits, so that the corep modes are combination 
of irrep modes on the two orbits. In this case, the two modes generated by IA3I are linearly 
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TABLE III: Magnetic (axial vector) corep modes for the perovskite i?-site (Mn in the case of 
REMnO^) associated with the four type-a irreducible co-representations for space group Pnma 
and propagation vector = (/i, 0, 0). Note the similarity of these modes with those listed in refj^i, 
(Tab. Ill), and references cited therein. However, through corep analysis, we have specifically 
enforced invariance by application of the antiunitary operator KI, where / is the inversion at the 
origin of the coordinate system and K is the complex conjugation. This invariance can only be 
accidental in irrep analysis. The matrix elements for unitary and antiunitary operators can be 
found in Tab. [H e = e"''*'^^ and e* is its complex conjugate. 







Mn(l) = i 0,0 
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independent. Tab. [IV] lists the corep modes obtained by symmetrizing the conventional 
irrep modes, following the procedure described in Eq. IA3I 

B. BiMn205 and DyMn205 

The space group is Pbam (no. 55 in the International Tables^), with three relevant 
Wyckoff sites: 4f (0, |,z) with z «i | for the Mn^~^ sites, 4h {x,y,^) for the Mn^^ and 
4g (x, y, 0) for the RE sites. 4h and 4g have the same symmetry, and can be treated in a 
completely analogous way. The propagation vector is (|,0,0), in Kovalev notation, for 
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TABLE IV: Magnetic (axial vector) corep modes for the perovskite ^-site {RE in the case of 
REMnO^) associated with the four type-a irreducible co-representations for space group Pnma 
and propagation vector kj = (/i, 0,0). By construction, these modes are invariant by application 
of the antiunitary operator KI, where / is the inversion at the origin of the coordinate system and 
K is the complex conjugation. The matrix elements for unitary and antiunitary operators can be 
found in Tab. [He = e"''*'^^ and e* is its complex conjugate. 





RE{1) 


= xA,z RE(2) 




-z + i RE{2,) = -x + i,|,z-i REU) = 


-X + 1, -4 




y 
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the low-temperature phase of Dy^^ (the minority component (|,0,/i) will be dealt with in the 
next section) and (|,0,|), /C24 for Bi"^^. Both propagation vectors are special points of the 
Brillouin zone for which k = —k, and have identical irrep/coreps. Moreover, as explained 
in Section [Tll the inversion operator / G G^, and there is no need to introduce antiunitary 
operators. The analysis is therefore very similar to the one performed by Munoz^^ for 
the Bi case. It is convenient to perform a unitary transformation of the Kovalev matrix 
representatives and associated modes, through the unitary matrix 



2 



1 i 



With this transformation, the resulting matrices (Tab. |V|) become real. For site 4f, the 
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modes (Tab. IVl|) can also be made real by multiplication of each subspace basis set by an 
appropriate coefficient. For site 4h, a bit more care is required, since each of the two irreps 
appears twice for every spin direction, so there is a degree of arbitrariness in the definition 
of the invariant subspaces. Here, we have chosen the definition so that the basis vectors 
have constant moments on all sites, but other choices are possible. Physically, this means 
that the magnetic moments are related in pairs (SI is related to S2 and S3 to S4), but the 
pairs are allowed in principle to have different moments within the same irrep. 

TABLE V: Matrix representatives of the irreducible representations of the little group G'^ for the 
space group G = Pbam and k2o = (|, 0, 0), or = (|,0,|). The matrices reported herein are the 
same as in Ref.^^, and are related to the Kovalev matrices by the unitary transformation UMU~^, 
where U is given in Eq. [71 



Irreps 



hi h2 hs h4 h25 h26 h27 h 



28 




The key aspect to assess the symmetry of the possible magnetic structures is the fact 
that both irreps are 2- dimensional. Therefore, in contrast with the REMnOs example, 
it is possible to obtain non-centrosymmetric and polar structures even for a single order 
parameter, provided that certain special directions in the 2-dimensional space are avoided. 
The magnetic basis vectors for the propagation vector A;24 = (|,0,|) are reported in Tab. IVII 
and Tab. IVIII for the two sites. Modes for = (|, 0, 0) are slightly different using our atom 
conventions, but can be obtained in the same straightforward way (see Table captions). 

The experimentally determined magnetic structure belongs to the Fi irrep, and is a 
combination of and my modes for the 4f sites and ml^ and ■m2y for the 4h sites. The 
important point is that only one unprimed basis vectors of each subspace is employed in 
each instance, so that all the components transform in the same way: 
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TABLE VI: Magnetic (axial vector) irrep modes for the 4f (0, \,z) sites {Mn^+ in REMn205) of 
space group G = Pham and /c24 = The basis vectors reported herein are related to the 

Kovalev basis vectors (obtained by the projection method with the help of the program Sarab2I by 
the unitary transformation vU~^, where U is given in Eq. [71 The matrix elements can be found 
in Tab. |Vl Primed and unprimed modes (e.g., rux and belong to the same invariant subspace, 
and have been both multiplied by the same coefficient, so that the unprimed mode is always 1 on 
atom 1. 
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ip{Mn'^^) = CiTTix + C2m'y = Ci [m^, ■ 

ip{Mn'^^) = c^mlx + C4m2y = C3 [ml2,,ml^] ■ 

The matrices are applied precisely to the column vectors in Eq. [HI With this notation, 
the symmetry can be read straightforwardly from the irrep matrices in Tab. [Vl It is clear 
that only diagonal matrices (/13, h2e, h2s) survive, because off-diagonal matrices transform 
unprimed into primed modes. The structural point-group symmetry in the magnetically 
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+ C2 [my, m'y] 



1 




(8) 
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+ C4 [m2y,m2'y] 



1 
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TABLE VII: Magnetic (axial vector) irrep modes for the 4h {x, y, ^) Mn^+ sites (x ~ 0.1, y ~ 0.85) 
of space group G = Pbam and /c24 = The RE atoms are on site 4g with the same site 

symmetry, and their modes can be deduced in a completely analogous way. The basis vectors 
reported herein are related to the Kovalev basis vectors (obtained by the projection method with 
the help of the program Sarah^I by the unitary transformation vU~^ , where U is given in Eq. [71 
The Sarah modes were preliminary recombined across invariant subspaces so as to have constant 
moments on all sites. The matrix elements can be found in Tab. |Vl Primed and unprimed modes 
(e.g., mix and ml'^) belong to the same invariant subspace, and have been both multiplied by the 
same coefficient, so that the unprimed mode is always 1 on atom 1. 
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1 



ordered phase is therefore m2m, which allow the polarization only along the b axis, as found 
experimentally. 

C. TbMn205 - commensurate phase 

TbMn205 is isostructural to the previous compounds, but orders magnetically with 
different propagation vectora^"^. Here, we will only consider the commensurate, high- 
temperature phase, with (|,0,|), kiQ in Kovalev notation, but the same analysis would 
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apply to an incommensurate phase of the type (|,0,/i), also labeled kiQ. There are only 3 
elements in the irrep little group: a two-fold axis || to z (/14) and two glide planes _L x (/126) 
and y (/?.27)- There is a single 2-dimensional irrep. The (3 matrix is the identity, so the 
matrices for the antiunitary operators Kh25, Kh2, Khs and Kh28 are the same as those for 
hi (the identity), h2e, /i27 and h^, respectively. The basis vectors can be constructed using 
Eq. IA3t and are obviously invariant by inversion. The 4f (0, i, \) (Mn^+) sites are split by 
the irreps into two orbits, which are recombined to obtain the basis functions for the coreps, 
whereas site 4h {x,y, |) remains as a single orbit and Eq. IA3I produces sets of degenerate 
vectors. At this stage, it is useful to perform a unitary transformation using the matrix 

which makes all the matrices real. Note the special form of the unitary transformations 
for antiunitary operators (Eq. IA2|) . The resulting corep matrices are reported in Tab. IVIIIi 

TABLE VIII: Matrix representatives of the irreducible co-representation for the space group G = 
Pbam and kiQ = (^, 0, /i) (/U = ^ for the commensurate phase of ThMn20^). The Kovalev matrices 
(the same for pairs of unitary and antiunitary operators, as explained in the text) were transformed 
using the unitary matrix from Eq. [9l 



Coreps 



h4 h26 h27 Kh25 Kh2 Kha Kh; 
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1 
1 



The basis vectors can be derived in the same way, by applying the inverse unitary trans- 
formation to the corep basis from Eq. IA3t and are reported in Tab. IIXI and 1X1 

It is easy to show that the constant-moments experimental solution^^i^Ei^ involves mixing 
modes spanning the same subspace with equal or opposite coefficients. Only the mix/ml'^ 
and to a lesser extent the mly/mVy components are relevant on both sites, so that, for 
example, the x-components ipx of sites Mn^+ and Mn^+ are: 
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TABLE IX: Magnetic (axial vector) corep modes for the 4f (0, l,z) sites (Mn^+ in REMn20^) 
of space group G = Pbam and A;i6 = (^,0,|). Primed and unprimed modes (e.g., mix and ml^ 
belong to the same invariant subspace. 
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'J* {mix + ml'x) = 7* [ml^.,ml^] ■ 



(10) 



V'x(Mn^+) = (ml, + ml^) = [ml„mi;.] 



where 7 



V2 



Once again, the matrices are apphed precisely to the column vectors in 



Eq. [TOl With this notation, the symmetry can be read off directly from Tab I Villi Note that 
the phase factor 7* in Eq. [10] does not reduce the symmetry even in the commensurate case, 
because 7^ is a lattice translation. The surviving operators h26, Kh3 and Kh28 define the 
structural point group symmetry m2m {C2v), indicating that b is the only polar direction, as 
observed experimentally. Note that the conclusion is identical to the BiMn205 case, in spite 
of the fact that we have adopted different basis conventions for the invariant subspaces. This 
is a direct confirmation of our conjectured^ that the c-axis component of the propagation 
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TABLE X: Magnetic (axial vector) correp modes for the 4h (x, y, ^) Mn^+ sites (x ~ 0.1, y ~ 0.85) 
of space group G = Pbam and kiQ = (^,0,|). The RE atoms are on site 4g with the same site 

symmetry, and their modes can be deduced in a completely analogous way. Primed and unprimed 

1 + i 

modes (e.g., mix and ml^.) belong to the same invariant subspace. 7 = — and 7* is its complex 
conjugate. 
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vector does not affect the symmetry properties of the system, provided that the in-plane 
magnetic structure is the same. 

D. HgCr2S4 

Recently, attention has been drawn to chalcogenide chromium spinels (cubic, space group 
FdZm) of the type ACr2X4 (A=Cd,Hg; X=S,Se) which are weakly ferroelectric in their 
magnetically ordered state and have been classify as multiferroic materials^i^i^. HgCr2S4 
is particularly interesting, because it has a complex magnetic structure, whereas all the 
other chalcogenide spinels are ferromagnetic. Very recently^, we have studied the HgCr2S4 
magnetic structure using high-resolution neutron powder diffraction. Long-range incom- 
mensurate magnetic order sets in at Tjv ~22K with propagation vector k=(0,0,~0.18). On 
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cooling below T^r, the propagation vector increases and saturates at the commensurate 
value k=(0,0,0.25). The magnetic structure below Tat consists of ferromagnetic layers in 
the a6-plane stacked in a helical arrangement along the c-axis. We also performed a full 
symmetry analysis using co-representations, determining the matrices and modes for the 
relevant corep, which is derived from the irrep as explained above. We will not repeat 
the detailed analysis, referring instead to our previous paper—, of which we summarize the 
salient point herein. There are 4 corep basis vectors, corresponding to pairs of ferromagnetic 
and antiferromagnetic helices, all with the moments in the x — y plane. Pairs of basis vectors 
defining invariant subspaces correspond to left- and right-handed helices. The experimental 
solution is one of the ferromagnetic helices: as in the previous case, the symmetry for an in- 
commensurate propagation vector can be read off the matrices in Ref.— , Tab. 11, considering 
that only those in diagonal form {hi, h^, hi^, /lis, Kh2, Kh^jKhi^ and Khi^) survive below 
the ordering temperature. These operators define the structural point group 422, which is 
non-polar, in spite of being non-centrosymmetric. We therefore conclude that ferroelectric- 
ity cannot arise from the magnetic transition at Tat ~22K to an incommensurate phase, 
although the observed polarization can very well have other causes. The low-temperature 
phase is close to be perfectly commensurate with k=(0,0,0.25), and the assessment of its 
symmetry requires additional care, since, as already illustrated, it depends on the overall 
phase factor of the magnetic structure. As a preliminary observation, we remark that the 
lattice is F-centered, so lattice translations have eigenvalues that are multiple of e**. The 
4-fold rotations hu and hi^ are always lost, because ie = e^s'^ is not a lattice translation. 
The symmetry is always lowest for a generic phase factor, for which only survives (point 
group 2). However, pairs of orthogonal 2- fold axes survive for a global phase of mod | 
{Khi3 s-iid Khio) or mod | {Kh2 and Kh^), so that the point-group symmetry is 222 in 
both cases. 



IV. DISCUSSION AND CONCLUSIONS 



We have presented a general method, based on a non-conventional application of co- 
representation analysis, to determine the point-group symmetry below a magnetic ordering 
transition for a crystal that is centrosymmetric in the paramagnetic phase, regardless of the 
direction and magnitude of the magnetic propagation vector. This method was employed 
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to determine the constraints on the development of electrical polarization in a new class of 
"magnetic improper" multiferroics. This approach can be readily extended, with essentially 
no modifications, to paramagnetic crystals with non-polar, non-centrosymmetric groups. It 
is also noteworthy that the point group we derive can be employed to set constraints on other 
vector or tensor quantities in the magnetically ordered phase. Pragmatically, we found the 
process of deriving corep matrices to be straightforward, when one has become familiar with 
the Kovalev tables^-. The only caveat is that one should be careful in dealing with Eq. IA41 
since the operator aQga^"^ may be related to by a translation, entailing an additional phase 
factor in the matrices. We found that the most effective way of deriving corep basis vectors 
is to transform the irrep basis vectors as in Eq. IA5I The irrep basis vectors can be obtained 
directly by projection, or with the assistance of one of several dedicated programs such as 
FullProF, Mody or Sarab2Ii2^i2£. Whatever the method employed, it is useful to check the 
symmetry of the basis vectors against the corep matrices. This is best done graphically, by 
superimposing the basis vector pattern, with the appropriate phase factors indicated, onto 
a diagram of the symmetry elements as in the International Tables^. With a savvy choice 
of the basis vectors, the structural point-group symmetry of the magnetically ordered phase 
can often be deduced by inspection for a variety of physically relevant magnetic structures^. 
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APPENDIX A: AN OVERVIEW OF CO-REPRESENTATION THEORY 

Coreps are constructed in a very similar way to irreps, on the basis of the subset of 
{E, K} G that leaves k invariant, the equivalent of the "little group" G^. It is important 
to stress that, by its very construction, (like G^) never mixes the +k and — k Fourier 
components, so the basis functions obtained through corep analysis are trivial eigenvalues of 
the pure translations. We can therefore drop the prefix e*'' * in Eq. [3] and limit our analysis to 
the representative elements of as for usual irreps. This accomplishes our goal of capturing 
the full symmetry properties of m whilst conveniently working with Fourier components. We 
can distinguish three cases: 

1. = G^. This occurs when — k is not in the irrep star of k, which is the case only 
for certain non-centrosymmetric space groups. By their very definition, the magnetic 
improper multiferroic materials we deal with here are always centrosymmetric in the 
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paramagnetic phase, and this case is therefore never relevant. 

2. = G^ + KG^. Since Kk = -k and KE e M^, k and -k must be equivalent, i.e., 
k must be a Lifshits vector. The use of coreps resolves the difficulty with inequivalent 
complex conjugate representations, which in the conventional irrep analysis are some- 
what artificially combined in "physical" irreps. However, coreps do not add much in 
terms of symmetry analysis, since all the relevant proper and improper rotations are 
already contained in G^. 

3. = G^ + aoG^, where = Kh and h E {G — G^). Once again, the very nature 
of our problem dictates that the inversion / G G, so we can always choose = KI. 
In this case, it is clear that contains more rotations than G^, and coreps should 
always be used. 

The small coreps of are built out of pairs of small irreps of G^, A{g) and A{g) = 
A(aQ^gaQ)* and their corresponding basis vectors ip and (j) = aQ-ip. Note that ijj and (f) may 
be linearly dependent, but are always independent if the magnetic atom is split into orbits 
or if A{g) and A{g) are not equivalent. Using the compound basis the derived corep 

is always of diagonal form for g G G^ and of off-diagonal form for a^g G {M^ — G^) : 



The matrices in Eq. lAll may be reducible to a simpler form by a change of basis through 
a unitary operator U. Note that the similarity transformation has a different form for 
antiunitary operators: 



A{g) 
A{g) 



(Al) 



D{aog) 



A{aogao) 
A{gy 



D'ig) =UD{g)U-' \/g E G'' 
D'iaog) = U*D{aog)U-' \/aog G (Af^ - G"^) 



(A2) 



We can distinguish three further cases: 
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A{g) and A{g) are equivalent through a unitary matrix A^, so that A{g) = NA{g)N~^ 
and NN* = +A(aQ). In this case, Eq. lAll can be reduced to block diagonal form 
for both unitary and antiunitary operators. A significant simplification occurs if A^ is 
the identity matrix. This is always the case, for example, when there is at least one 
g G that commutes with and has a real matrix representative. With the new 
basis system: 

^' = ^^ + 0, ij" = i{ij-(P) (A3) 
the corep decomposes into two identical coreps: 

gij' = A{g)ij\ aogf = A{a^ga^^)^' (A4) 
gr = A{gW\ aogr = A{aoga^')r 

The second corep is often made antisymmetric with respect to the antiunitary op- 
erators by omitting the imaginary unity in the construction of ip"- Clearly, the two 
resulting coreps and D_ remain identical. If A^ is not the identity, Eq. IA3I must 
be generalized to: 

ijj' = ^Ij + N*(f), ^" = i{'4) - N*(l)) (A5) 
All the examples in Section IIIII belong to this " Case a" . 

A{g) = NA{g)N~^ but NN* = — A(aQ). In this case, which is comparatively rare for 
centrosymmetric groups, Eq. lAll cannot be reduced to a diagonal form. Instead, an 
appropriate transformation is applied to convert the matrix representatives of antiu- 
nitary operators into a block-antisymmetric matrix. 

A{g) and A{g) are not equivalent. In this case, we retain the form of Eq. lAll with the 
same basis. Note that A{g) must necessarily be equivalent to one of the other irreps 
in the list, say A{g) = NA'{g)N-'^. 
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Complete tables of the coreps for all crystallographic space groups, as well as of the 
"auxiliary matrices" (therein called (3) are contained in Ref.— . From these tables, one 
can readily construct the corep matrices and the new basis vectors. 
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